Non-sequential double ionization of molecules in a strong laser field 
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We consider the final stage of double ionization of O2 molecules by short linearly polarized laser 
pulses. The saddles of the effective adiabatic potential energy close to which simultaneous escape 
of electrons from a molecule takes place are identified. The analysis of the saddles and numerical 
simulations of the ionization indicate that to observe clear signatures of simultaneous electron escape 
in double ionization of O2 molecules stronger and much shorter laser pulses than those used in the 
recent experiment [E. Eremina, el al, Phys. Rev. Lett. 92, 173001 (2004)] should be applied. 
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I. INTRODUCTION 

Experimental studies of a non-sequential double (or 
multiple) ionization is possible thanks the use of high 
intensity ultrashort-pulse lasers 0,0- On the other hand 
theoretical description of such multi-electron dynamics is 
far from complete and suggests that interactions between 
electrons become important in considerations Q, 0, IE H| • 
In this paper we present description of the non-sequential 
double ionization of molecules based on a classical model 
for electrons in a combined Coulomb and external field 
following the approach developed in 0, IE IE for the 
multiple ionization of atoms. 

Recent experimental investigations aimed on the non- 
sequential double ionization of diatomic molecules [Tol 
ITU Il2l Il3| showed that there are differences between 
molecular species. For instance, in the case of N2 it 
seems that electrons escape with similar momenta along 
field polarization axis more often than in the case of 
O2 |l3j |. Such a correlation in momenta of escaping 
electrons is viewed as characteristic feature of the non- 
sequential double ionization. In our previous publication 
|l 1| we have analyzed non-sequential double ionization 
of molecules and presented numerical results mostly for 
double ionization of N2 molecules. Here we complete the 
analysis by investigating how the distributions change in 
the case of O2 molecules for different parameters of the 
system. 



II. RESULTS 

In our considerations we assume the re-scattering sce- 
nario [T5I fl6| . Namely, in the first step one electron tun- 
nels out through the Stark saddle and then is returned 
back to the nucleus. And so at the expense of the energy 
brought back by the returning electron a highly excited 
state of a molecule is formed. Finally, such a highly ex- 
cited compound state decays in several ways through a 
single, double or multiple ionization. We start our anal- 
ysis after the formation of the excited complex that is to 
say we assume that we have an initial state of two highly 



excited electrons close to the molecular core in the pres- 
ence of linearly polarized laser field. For short laser pulses 
molecules have not enough time to change their orienta- 
tion 0, 0| and thus we presume that the motion of the 
molecular core is frozen. Hence, the Hamiltonian reads 
(in atomic units, which are used throughout): 
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The potential consists of the potential energies associ- 
ated with interactions of the electrons with the nuclei 
(the entire structure of the molecular core is approxi- 
mated by two positively charge nuclei), with each other 
and with the external field. Origin of coordinate system 
is placed in the center of mass of the nuclei and without 
loss of generality it is assumed that the molecule lies in 
the x-z-plane. Then n, yi and Zi indicate position of the 
electrons, d, the distance between nuclei, and 6, the angle 
between the molecular axis and the z-axis (the polariza- 
tion axis), respectively. The electric field strength F(t) 
has an oscillatory component times the envelope from the 
pulse, namely: 



F(t) = Fcos(ujt + 0) sin 2 (7rf/T d ), 



(3) 



with F, w, <j) and Td the peak amplitude, frequency, initial 
phase of the field and pulse duration, respectively, 
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where n is number of cycles in the pulse. 
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Within this model the only difference between differ- 
ent diatomic species lies in the distance, d, between the 
nuclei. In the present publication we restrict ourselves to 
the oxygen molecule, hence d = 2.28 a.u. in the following. 

To identify channels for simultaneous electron escape 
we use an adiabatic approximation, keeping the field 
fixed, for the reason that the classical motion of the elec- 
trons is fast compared to the field oscillations (later, in 
the numerical simulations, we use the potential @ with- 
out any adiabatic assumption) . The channels for the non- 
sequential ionization correspond to saddles of the poten- 
tial (J2J- That is, when electrons pass close to a saddle 
they may leave a molecule simultaneously because one of 
unstable directions of a saddle allows both electrons to 
increase distance from a molecule. 
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FIG. 1: Energy of the saddle (top panel) and cross section 
exponent (bottom panel) as a function of orientation angle, 
8, for O2 molecule (d = 2.28 a.u.). Broken line corresponds to 
the saddle in the yz-plane, whereas solid line corresponds to 
the saddle in the :rz-plane. The field strength is F(t) = 0.075 



For a general orientation of a molecule with respect to 
the field polarization axis we have found two saddles for 
electron escape. The first is situated in a plane defined 
by the molecular and field axes while the other in the per- 
pendicular plane. A few stable and several unstable di- 
rections are revealed by the local stability analysis of the 
saddles. One unstable direction corresponds to the non- 
sequential double ionization, while the other unstable di- 



rections reflect interactions that will push electrons away 
from the non-sequential double ionization path leading 
to the single ionization. 

The Lyapunov exponents that characterize the differ- 
ent directions allow one to obtain the cross section behav- 
ior for non-sequential double ionization events. Similarly 
to the double ionization without a field, analyzed many 
years ago by Wannier |18t Il9l 12(1 l2l| , the competition 
between the various unstable directions gives rise to an 
algebraic variation of the cross section with energy close 
to the threshold, namely, 



a(E) <x(E- Vs) c 
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where Vs is the saddle energy and the exponent contains 
the Lyapunov exponents, 



(6) 



A r is the Lyapunov exponent of the unstable direction 
corresponding to the non-sequential double ionization 
path, and are the Lyapunov exponents of all other 
unstable directions of the saddle 0, E3 • 

In Fig. H variations of the parameters of the saddles 
with the orientation 9 for the internuclear distance corre- 
sponding to O2 molecule are presented. Fig-Hshows that 
the energies and exponents do not change very much with 
the orientation of the molecule (especially in the case of 
the saddle located in a plane defined by the field and 
molecular axes). Thus, we may conclude that from the 
point of view of the local analysis of the non-sequential 
ionization, decay of the highly excited complex does not 
depend strongly on the orientation of the molecule with 
respect to the field polarization axis. 

Thus far we have analyzed the potential within the adi- 
abatic approximation. Now we present results of numer- 
ical simulations of decay of a highly excited two electron 
compound state taking into account full time dependence 
of the potential . Before the double ionization escape 
both electrons pass close to the nuclei where they interact 
strongly with each other and with the nuclei. Thus we 
may assume that all memory of the previous motion (i.e. 
tunneling of the first electron, its evolution in the com- 
bined Coulomb and laser fields and the rescattering pro- 
cess) is lost. Then it is plausible to assume that the com- 
pound state which decays to a doubly charged molecule 
can be classically modeled by a statistical distribution for 
two electrons close to the nuclei. Detailed discussion on 
numerical calculations and initial conditions is presented 
in our previous publication fl4j |. 

In Fig. [5] distributions of the final electron momenta 
parallel to the field axis for the O2 molecule oriented 
along the field axis are presented. Different initial en- 
ergy E and different pulse durations are considered. For 
very short laser pulses (n = 2 cycles) the distributions 
are localized along the diagonals indicating that the elec- 
trons escape predominately by passing close to the sad- 
dles analyzed previously. For longer pulses and for energy 
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FIG. 2: Final distribution of the parallel electron momenta for different initial energy, E and different pulse length (i.e. number 
of cycles in the pulse, n) in the double ionization of the O2 molecule. The field strength is F — 0.075 a.u., the frequency 
Lu — 0.057 a.u., the internuclear distance d = 2.28 a.u. and the orientation parallel to the field, — 0. The plots differ in initial 
energy E and the number of field cycles n: (a) E — —0.05 a.u., n = 2; (b) E — —0.05 a.u., n = 26; (c) E = —0.9 a.u., n = 2 
and (d) E = -0.9 a.u., n = 26. 



E = —0.9 a.u. the distribution drastically changes its 
character. Namely, the first and fourth quadrants of the 
panel become strongly populated implying that a number 
of sequential decays significantly increases. The reason 
for this is that, after re-scattering when a highly excited 
two electron complex is created, there are two dominant 
scenarios within the first half cycle of the field: a non- 
sequential double escape or a single ionization. During 
the next cycles, unless the molecule is already doubly 
ionized, we are left with a singly ionized molecule which 
may survive to the end of the pulse or the second electron 
escapes and that corresponds to a sequential double ion- 
ization. As the pulse becomes longer, the number of such 
outcomes increases, even to the point of overwhelming 
the number of direct, non-sequential double ionization 
cases. Hence, even though the re-scattering scenario is 
involved in the double ionization process, for sufficiently 
long pulses the momentum distributions will not show 
signatures of the non-sequential electron escape. For the 
initial energy E much higher than the minimal energy 



of a saddle (Vs ~ —1.27 a.u.) the probability of non- 
sequential ionization is bigger and even if the pulse du- 
ration is quite long (e.g. n = 26 in panel (b) of Fig. [2J 
the signatures of non-sequential process remain. 

Alteration in the orientation of the molecule does not 
lead to significant changes in the momentum distribu- 
tions as one can see by comparing Fig. [3] and Fig. For 
the field amplitude used in the simulations the saddles for 
non-sequential process are far away from the nuclei and 
the positions and other parameters of the saddles change 
only slightly with the vary of orientation 9 — compare 

Figs. 12131 should be compared with the corresponding 
figures for nitrogen molecule presented in Ref. [l4| - 
there we chose energies E = —0.3 a.u. and E = —0.6 a.u. 
(also the field peak amplitude was slightly smaller, i.e. 
F = 0.07 a.u.). Comparing the data we may conclude 
that high initial energy leads to much clearer signatures 
of simultaneous electron escape (even for long pulses) 
than low energy. Initial energy is determined in the 
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FIG. 3: Final distribution of the parallel electron momenta 
for the O2 molecule at non-zero angle to the field axis, i.e. 
9 — 7r/4, and for different initial energy E. F — 0.075 a.u., 
ijj — 0.057 a.u., d = 2.28 a.u. and n = 2. The energies are (a) 
E = -0.05 a.u. and (b) E = -0.9 a.u. 



rescattering process, thus to enhance signatures of the 
simultaneous escape one should apply stronger external 
field. Increasing field amplitude is, however, dangerous 
because for very strong laser pulses sequential ionization 
(without rescattering) becomes dominant. Compromise 



can be reached by increase of the amplitude with simulta- 
neous reduction of the pulse duration — very short laser 
pulses and the field stronger than that used in Ref. 
should allow for clear observation of simultaneous elec- 
tron escape in double ionization of O2 molecules. 

Our analysis is restricted to the final stage of the non- 
sequential double ionization of O2 molecules, where we 
choose initial conditions microcanonically. Early stages 
of the excitation process and the nature of the compound 
state (i.e. the symmetry of the valence orbital may af- 
fect the excitation process) before the final decay towards 
double ionization may have influence on the initial con- 
ditions. That, in turn, may modify the final parallel mo- 
menta distributions. This aspect needs further investiga- 
tion. 



III. SUMMARY 

We have performed a purely classical analysis of the 
final stage of the non-sequential double ionization of O2 
molecules in the strong laser field where a molecular core 
is approximated by two positively charged nuclei and the 
re-scattering scenario is assumed. Our analysis, together 
with the previous studies devoted mostly to double ion- 
ization of N2 molecules [Ti}. allow us to draw following 
conclusions: i) From the point of view of classical anal- 
ysis, within the considered model, there is no difference 
between nitrogen and oxygen molecules in the sense that 
both of them can show signatures of simultaneous dou- 
ble escape, ii) Orientation of the molecule with respect 
to the field axis does not influence significantly the fi- 
nal momentum distribution for the initial energy range 
considered in the model, iii) Shorter and stronger laser 
pulses should enhance signature of simultaneous electron 
escape that will be visible in the momentum distribu- 
tions as a more pronounced symmetrical escape of the 
electrons. 
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